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We analyze the dynamics and the geometric properties of the Potential Energy Surfaces (PES) 
of the fc-Trigonometric Model (fcTM), defined by a fully-connected fc-body interaction. This model 
has no thermodynamic transition for fc = 1, a second order one for fc=2, and a first order one 
for k>2. In this paper we i) show that the single particle dynamics can be traced back to an 
effective dynamical system (with only one degree of freedom); ii) compute the diffusion constant 
analytically; Hi) determine analytically several properties of the self correlation functions apart from 
the relaxation times which we calculate numerically; iv) relate the collective correlation functions 
to the ones of the effective degree of freedom using an exact Dyson- like equation; v) using two 
analytical methods, calculate the saddles of the PES that are visited by the system evolving at 
fixed temperature. On the one hand we minimize |Vyp, as usually done in the numerical study 
of supercooled liquids and, on the other hand, we compute the saddles with minimum distance (in 
configuration space) from initial equilibrium configurations. We find the same result from the two 
calculations and we speculate that the coincidence might go beyond the specific model investigated 
here. 

I. INTRODUCTION 

In the last few years there has been an intensive study Q, Q, 1^ ^ IE IE Q concerning the connection between 
the slow dynamics of complex systems and the extrema of the potential energy surface (PES), defined through 
the potential energy function V{q), following the seminal work of Stillinger and Weber 

Two ways of studying the dynamics of supercooled liquids and glasses that are based on the analysis 
of the PES have been of particular importance in the recent past. The first one concerns the analysis of 
the properties (energy location, number, curvature,..) of the minima of the PES that are visited by the 
system during its evolution at fixed temperature. Assigning to each minimum its zero-temperature basin of 
attraction, one obtains a partition of the phase space allowing to define a configurational entropy for the 
supercooled and the out-of-equilibrium glassy regime . The properties of the minima of the PES have also 
been connected to several features of supercooled liquids and glasses. We can mention their relation to the 
fragility of the glass former 3] , the diffusion processes in supercooled liquids 0, 0, ll^- and the effective 
fluctuation-dissipation temperature Iislll^ in the out-of-equilibrium glassy phase |l^. This method 
is closely related to Edwards' proposal [1^ to describe the main properties of granular matter with a flat 
measure over blocked configurations that correspond to the minima of the PES. Note that granular matter 
is an effectively zero-temperature system for which the study of the energy landscape if fully justified. 

The second approach, which corresponds to the study of the instantaneous normal modes, is based on 
the study of the eigenvalues of the Hessian at the instantaneous configurations along the trajectory (in 
configurational space) that the system follows during its dynamical evolution 0,153 ■ This approach allowed 
one to relate the diffusion process to the local curvature of the landscape. 

More recently a third approach has been proposed li^,^ and apphed j^E, IH, ^5*, '26*, '2^, HI] to study the 
slow dynamics in supercooled liquids. Within this approach, the saddles of the potential energy surface play 
a central role. It has been found numerically that the order of the saddles (number of negative eigenvalues 
of the Hessian matrix) visited during the equilibrium dynamics at temperature T extrapolates to zero when 
T reaches the dynamic transition temperature Tmct (or mode-coupling temperature (29j ) . 
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The role of the stationary points of the PES (saddles and minima) has been also pointed out in a different 
context. Indeed, studies aiming to clarifying the microscopic origin of phase transitions suggest that the 
presence and order of such transitions is related to changes in the topology of the manifold of the PES 
sampled by the system when crossing the (thermodynamic) critical point [sol l3lj | . This has been observed 
by counting the number and the order of the stationary points of V{q) and building up the Euler characteristic 
of the manifold. The latter is a genuine topological property of the energy surface defined by V{q) ^constant, 
and, in particular, it does not depend on the statistical measure defined on it (i.e., on temperature). 

Disordered mean-field spin models have been proposed to mimic the behaviour of super-cooled liquids and 
glasses. Their statics and dynamics, as well as the properties of their free-energy and energy landscapes, are 
amenable to analytical studies p^ . The main features mentioned in the previous paragraphs are realized 
by these models where, at finite temperature, the geometry of the free-energy landscape replaces the PES. 
In particular, the importance of saddles in the free-energy landscape for the evolution of these systems has 
been elucidated in the past and a comparison between the roles played by free-energy and energy landscapes 
has also been discussed in this and more general contexts 32, 33]. 

If one wishes to relate the equilibrium dynamics of a complex system to the properties of the saddles of its 
PES, an unambiguous definition of saddle visited during the equilibrium dynamics is mandatory. Until now, 
two different definitions have been used: 1) In the numerical simulations of simple models - but still too hard 
to study analytically - such as Lennard- Jones systems, a partitioning of the configuration space in basins 
of attraction of saddles is obtained via an appropriate function W (usually W = |Vyp) that has a local 
minimum on each stationary point of and the saddles are then obtained via a minimization of W starting 
from an equilibrium configuration obtained from a molecular dynamics simulation at temperature T. 2) In 
the analytic computations applied to disordered mean-field spin models one looks to the saddles that are 
closest, with respect to some distance in the configuration space, to a reference configuration extracted from 
the Gibbs distribution at temperature T 0| . Unfortunately, until now the two methods have been applied 
to different models so the comparison between them is still qualitative. 

In this paper we study a very simple mean-field model without disorder, the k Trigonometric Model 
(fcTM), for which one can calculate analytically all the relevant quantities that have been previously studied 
numerically for more realistic models. In spite of its simplicity, the thermodynamic behavior of this model 
is quite rich, and its PES shows some of the features that have been observed in the PES of Lennard- Jones 
systems 0, [23, H^, 13 US HE 113 ■ Unfortunately, the model is too simple to show an interesting dynamics. 
The dynamical behaviour is closely related to the thermodynamics and, due to the absence of frustration 
or disorder, no dynamical arrest is observed. Still, on the one hand we have been able to check analytically 
the validity of some ideas that had been proposed in the literature and, on the other hand, to elaborate 
a method for the minimization oi W = |VFp that will be of use for a larger class of mean-field models 
including the disordered ones. This method allowed us to compare analytically the two definitions of closest 
saddles to equilibrium configurations, and to show that they coincide in our model. 

The paper is organized as follows: in section^ we introduce the model and present its main features. 
Then we analyze its static properties: in section Hill we discuss its thermodynamical behavior, in section Hvl 
we study the topological properties of the PES and relate them to the results of section IIIII In sections Ivl 
and I VII we present a detailed study of the dynamical behavior of the model. Then, in section FVIII we discuss 
the definition of the closest saddles to equilibrium configurations, their properties and their relation with 
the dynamics of the system. Many of the calculations require the introduction of a formalism that may not 
be familiar to all the readers and is not really necessary to follow the relevant part of the presentation; they 
are then presented in detail in the appendices. 

II. THE MODEL 

The fc-Trigonometric Model (fcTM) has been introduced in with the aim to study the relation between 
phase transitions and topological properties of the PES. The model is defined by the Hamiltonian 



"=J^ E [1 - cos((p., + • • • + J] = iVA(l - SR^'^) , 



(1) 
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having introduced the "magnetization" 

i 

Here </JiG[0, 27r), ■ ■ ■ , N, are angular variables and A is the energy scale. It is easy to see that the model 
is not invariant under continuous transformations of ipi but only under the discrete group Ckv generated by 

^^^^' + T' (3) 

If one interprets the variable Lpi as the angle between a two-dimensional unitary vector and a fixed axis, the 
trasformations in Eq. (jSJ are rotations by an angle 27r/fc of the vector and the reflection with respect to 
the fixed axis. In the low temperature phase this symmetry is broken, and a spontaneous magnetization is 
generated in a direction %pn — 2TTn/k. We will often choose ip — without loss of generality, in order to keep 
unbroken the symmetry ipi ^ —ipi [s^. 

The system is subject to a dynamics generated by a Langevin equation j3^ : 

I'Pi = - J. 'r-q, , (4) 

where 7/ A is the time scale and rji is a Gaussian noise with 

{m{t)iij{t'))=2T-i5,j5{t^t') . 

We will consider averages of a generic observable yl({iy9}) over the noise distribution, and we choose random 
initial data (that corresponds to a quench from infinite temperature) . Eventually, we will consider the 
limit, in which the system equilibrates and is described by the Gibbs ensemble at temperature T. 



III. THERMODYNAMICS 



The termodynamics of mean field models is exactly solved neglecting the correlations between different 
degrees of freedom and obtaining an effective Hamiltonian that contains a parameter to be determined self- 
consistently. For example, in the fully connected Ising model, with Hamiltonian H = —{2N)~^ J^ij SiSj, the 
substitution reads SiSj — > {si)sj + Si{sj) — {si){sj). Defining m = (si), one obtains the effective Hamiltonian 
Ti.{s) = ms + c(m) for a single degree of freedom s (c(m) is an irrelevant constant that depends only on m). 
The self-consistency equation is finally obtained calculating m = (s) on this effective Hamiltonian. In fact, 
one can show that this procedure is equivalent to the evaluation of the free energy at the saddle point in the 
N —^ 00 limit. 

The generalization of this procedure to the fcTM is obtained substituting in Eq. ^ the expression: 

e'^"! . . . 6*'^-'= ^ k e'^'i (e*'^'2) • • • (e^'^''=) - (fc - 1) (e^'^'O • • • (e*'^-'=) . (6) 

and introducing the mean (complex) "magnetization" C=(e*'''), that has to be determined self-consistently 
on the mean field effective Hamiltonian H. As we always chose — 0, ( is real and the effective Hamiltonian 
reads: 

W = A[l + (fc-l)C''-fcC''"^cos<^] . (7) 
The self consistency equation for ( turns out to be: 

^=<^°^^^" = /o(/?AfcC^-^) ' 
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FIG. 1: Mean magnetization as a function of temperature for k = 1,2,3. For A; = 3 the value of the canonical 
transition temperature To = 0.63 is indicated by a dotted vertical line, while the temperature Tc = 0.72 at which the 
magnetic solution disappears is marked by a white dot. The magnetic solution is metastable for To < T < Tc; the 
same happens to the C = solution for T <Tq. 



where /3=1/T, Io{a) = {2Tr)~^ J^"" dipe°"=°'"^ and Ii{a) = I^ia) are the modified Bessel functions of order 
and 1 respectively. For each (3 the Eq. ((SJ gives the thermodynamic value of the mean magnetization 
C(r). The value C = always solves Eq. but is a stable solution only at low /? (high temperature). 
As P is increased other solutions may appear, and one has to consider the one that minimizes the free 
energy as the stable one, while the other solutions can be either unstable or metastable. In Fig. ^ we 
report the function C(^) for k = 1,2,3. For k=l the curve is smooth, no phase transition occurs and the 
magnetization approaches zero at high temperature. For fc = 2 a second order phase transition takes place 
at Tc = A, separating a high temperature paramagnetic phase, where only the solution C = exists, and a 
low temperature ordered phase, where C = becomes a maximum of the free energy separating two minima 
with C 7^ corresponding to "0 = and -0 = tt. From this symmetry structure one sees that for k ~ 2 the 
model is in the universality class of the Ising model; the (scalar) order parameter is the real part of the 
magnetization, while the imaginary part is never different from zero. The critical exponents are then the 
classical mean- field exponents of the Ising model; in particular we have C ~ |T — Td^/^ close to Tc and, if 
a perturbation 6Ti — —hcosip is added, C ^ h}/'^ at T = Tc- For fc = 3 (and also for fc>3, not reported in 
the figure) the system undergoes a first order phase transition. At high T only the paramagnetic solution 
C = exists, but on lowering T two other solutions appear at T = Tc- the magnetic one with C > ^-iid 
degeneration k (i.e. k different possible values of '0) and the one corresponding to the maximum of the 
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FIG. 2: Potential energy e (full lines) and saddle energy (dashed lines, see section lVnil as a function of temperature 
for fc= 1, 2, 3. The notation is the same as in the previous figure. 

free energy separating the magnetic and paramagnetic solutions. The magnetic minimum becomes stable 
at T = Tq < Tc- In Fig. ^ for the k = 3 case, we report the evolution of the two minima as a function 
of temperature. The second minimum appears at Tc — 0.72 (white dot in the figure) and the transition 
temperature Tq = 0.63 is indicated by a dotted vertical line. In Fig. |21the temperature dependence of the 
potential energy e = (7i) = A(l — (full lines in the figures) is shown for k — 1,2, 3. The first order phase 
transition also manifests itself in a discontinuity of the potential energy as a function of temperature. 

IV. TOPOLOGICAL PROPERTIES OF THE ENERGY SURFACE 

In this section we will study the properties of the stationary points (saddles) of the Potential Energy Surface 
(PES) of the system, defined by the Hamiltonian We will now focus only on the topological properties 
of the saddles, while in section IVIII we will study the properties of the saddles sampled by the system 
equilibrated at temperature T. We will now follow the derivation in [sj, while in appendix^ we present 
a different derivation that will be useful in the following. The stationary points (p are defined by the 
condition dH{ip)=0, and their order v is defined as the number of negative eigenvalues of the Hessian matrix 
llij{<f') — {d'^H/d(pid(pj)\^. To determine the location of the stationary points we have to solve the system 



— = - A fc ^[iz''-^e''P^] = A fc e''"^ sin[(A: - l)'tP + ipj] = , Vj , 



(9) 
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where we have used the definition z=£^e^^ given in Eq. A first group of solutions arises for ^ = 0; from 
Eq. we have H{(f)=NA[l — ^'^ cos(fc-0)], and then the stationary points with ^{(p)=0 are all located at 
the energy e=H{<f)/N=A. We will now restrict ourselves to the region e 7^ A because, as we will see at 
the end, the quantities in which we are interested are singular when e = A. The presence of this singularity 
seems to be related, as we discussed elsewhere [^H, to the presence (and the order) of a phase transition. 
For e ^ A, Eq. © becomes 

sin[(fc - 1)^ + ip,] - , Vj , (10) 

and its solutions are 



Lpf- = [mjTT - (fc - l)V']mod 277 , (H) 

where mj-ejO, 1} and m = {rrij}. Therefore, beside the different possible values of iJj, each stationary point 
(^™ is characterized by the set m. To determine the unknown constant ^p we have to substitute Eq. Hll|) in 
the self-consistency equation 

z = ^e'^ = iV-iy e^'c^ = (-1)™^' . (12) 

If we introduce the quantity n(Cp) defined by 

n = 7V-iy m, , 1 - 2n = iV-^V (-1)™^ , (13) 

and we have from Eq. H12|l 

e=|l-2n|, (14) 

^ |2Wfc forn< 1/2 , 

|(2Z + l)7r/fc for n> 1/2, ^ ' 

where / gZ. Then the choice of the set m is not sufficient to specify the set because the constant -0 can 
assume some different values. This fact is connected with the symmetry structure of the potential energy 
surface (the different values of V'z generate the multiplets of stationary points). We have then obtained that 
all the stationary points of energy e^^A have the form 

Cpf'^ = [m^TT - (fc - l)V';]mod 271- • (16) 

The Hessian matrix is given by 

H,j = A fc ^\N-^{k - i);2fc-2e«(¥'.+v.) + J^^.^fe-ie*'^')] . (17) 

In the thermodynamic limit it becomes diagonal 

Hy = b^j A fc i,^-^ COS {i){k - 1) + 'p^) . (18) 

One cannot a 'priori neglect the contribution of the off-diagonal terms to the eigenvalues of H, but we have 
numerically checked that their contribution changes the sign of at most one eigenvalue out of iV. Neglecting 
the off-diagonal contributions, the eigenvalues \j of the Hessian calculated at the stationary point Cp are 
obtained substituting Eq. ifT^ in Eq. (fH^ : 

= (-1)'"^ A fc i^-^ , (19) 



so the stationary point order ^{(p), defined as the number of negative eigenvalues of the Hessian matrix, is 
simply the number of rnj = l in the set m associated with ip; we can identify the quantity n{(p) given by 
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FIG. 3: The logarithm a of the number of saddles and the saddle order n as a function of the energy level e for 
fe = l,2,3. 



Eq. H13|) with the fractional order v{lp)/N of tp. Then, from Eqs. (^), H14|l and (|15|l we get a relation between 
the fractional order n{(p) and the potential energy e{(p)=H{(p)/N at each stationary point (p. It reads: 



"(e) = ^ 



1 - sg: 



i/fc' 



(20) 



Moreover, the number of stationary points of given order v is simply the number of ways in which one can 
choose V times 1 among the {rrij}, see Eq. ()16|l . multiplied by a constant that does not depend on N and 
takes into account the degeneracy introduced by Eq. H15|l . Therefore: i) the fractional order n=v/N of the 
stationary points is a well defined monotonic function of their potential energy e, given by Eq. (j2()|l . and ii) 



the number of stationary points of a given order v is Ck (^) . We can define the quantity 



1 f N 
aie) = lim — log L,r / - 

' JV^oo N \Nn{e] 



-n{e) log7i(e) - (1 - n(e)) log(l - n(e)) , 



(21) 



that represents the "configurational entropy" of the saddles. In jSJ we have shown that this quantity is 
related to the Euler characteristic of the manifolds = {'p\H{ip) < Ne} and that its singular behavior 
around the point e = A is related to both the presence and the order of the phase transitions that occur for 
k > 2. In Fig. |2|the quantity a{e) and n{e) are reported for k = 1, 2, 3 for all values of e 7^ A: one can see 
that the presence of a phase transition for k>2 is signaled by a singularity in the first derivative of a{e). The 
order of the transition seems to be related to the sign of the second derivative of cr(e) around the transition 
point, that is negative for second order phase transitions and positive for first order ones. 
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V. SINGLE-PARTICLE DYNAMICS 



The single-particle dynamics can be studied by means of observables of the form 

^(M)-^E^(^') ■ (22) 

i 

Some interesting quantities are, for example, the diffusion constant D(T)^ defined as 

D{T) = ^Um ^ - ipm?) , (23) 

i 

(where obviously the angular variables have to be considered as variables defined on the whole real axis 
without introducing the periodic condition ipi + 2t: = tpi) and the self correlation function, defined as 

i^(t,i') = i-^ [(e*['^'(*)-'^'(*')l) - (e*'^'(*))(e-"^'(*')) . (24) 

i 

To study the single-particle dynamics, we introduce an effective dynamical system for a single degree of 
freedom ip^ defined by 



im = -^{^{t),t)+v{t) , (25) 



where again 77 is a (5-correlated Gaussian noise with variance 2^T and the averages are calculated on its 
distribution. Note that the derivative of the effective Hamiltonian 7i((p, t) is performed with respect to tp 
at fixed time t (the time-dependence of TC is encoded in C{t))- The last equation is, as in the static case, a 
self-consistency equation. In appendix El we show that the averages of observables of the type H22|l on the 
dynamics defined by Eq. Q can be calculated using this effective single- variable dynamical system: 

{AiW})) = {Ai^))n . (26) 

We will consider the limit t ^ 00 in which the system is in equilibrium, and C does not depend on t and is 
equal to its equilibrium value ({(3) given by Eq. (jSJ. In this limit the correlation function (|24ll depends only 
on the time difference, F{t,t') = F{t — t'), and the dynamical system (|25|l reduces to 

795(0 + fcAC(/3)'=-i sin ^{t) =r^{t) . (27) 

If C 7^ 0, we can define the reduced variables i = — t and fi{t) = fc^^t-i v{t)i and Eq. (|?7jl can be 
rewritten as 

^+sin(^(f) =7?(i) , (28) 
at 

with {fi{i)fi{0)) — = 2fS{i). The k dependence is then encoded in i and f. 

A. The diffusion constant 

The analytical expression for the diffusion constant of Eq. 128|) is found for example in [s^l and is given by: 

b{f) = hm - ^(0)p) = — |— . (29) 



0.0 0.2 0.4 0.6 0.8 1.0 

T/A 



FIG. 4: Logarithm of the diflusivity D as a function of the temperature T for k — 1,2,3. For fc = 3 the notation is 
the same as in Fig.0 The dashed hne is the value of D in the paramagnetic phase, D = r/7. 



Then, from the definition of D given by Eq. (|23|) and from Eq. (|26|) one gets 

D{T) = ^ D{T) = -/o(fc/?AC(T)'=-i)-2 . (30) 

7 7 

In Fig. 2]we report Z3 as a fmiction of temperature T in a semilogarithmic scale. At low temperature Eq. (|30() 
predicts an Arrhenius behavior: 

D{T « A) ^ ^![^^£!le^/32Afe _ (3^) 
7 

In the high temperature paramagnetic phase one has C = so that Eq. (|28() reduces to a free Brownian 
motion and the diffusion constant is simply given by D — T/7 l38| . It is worth to note that for k = 1 the 
paramagnetic phase does not exist and D = T/j is only the asymptotic limit of Eq. H3U|) for T ^ 00. 



B. The self correlation function 



As we already discussed, the model is not rotationally invariant, and when the Ckv symmetry is broken a 
spontaneous magnetization appears, which phase can be an integer multiple of 27r/fc. We will choose the 
phase to be zero in order not to break the tpi — > —ifi symmetry. In this situation it is expected that even 
and odd functions of tp have different behavior. Therefore, it is convenient to separate the contributions 
in F{t) given by the real and imaginary part of e"^'. Using the relations {cos{(pi(t)) sin{(pi{0)) — and 
{sm{(pi{t))) = 0, due to the unbroken symmetry (pi —fi^ we can define from Eq. H24() (setting t' = Q 
because of the time translation invariance) : 

F{t) = F,{t) + F,{t) , 

■^c(i) = Tiy^ }{cos ip.i{t) cos ipi{{))) - (cos (^i(t)) (cos 93.^(0))] , 

N Y (32) 
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FIG. 5: The correlation function of costp and sin^ calculated using Eq. 1281 for two different values of T. They are 
well fitted by an exponential form. 



As previously discussed, see appendix IBl or Eq. H2t)|) . the above correlations are equal to the ones calculated 
using the effective dynamical system (|25|) : at equilibrium they are given by 



We will use for T{t) the same notation used for F{t)^ see Eqs. (I32f) . 



(33) 



1. The reduced system 

To compute the correlation functions, it is useful to use again the reduced variables and Eq. H28|) . Unfor- 
tunately, we have not been able to derive an analytic expression for the correlations in the whole T range, 
but only in the high and low temperature limits. In these limits the correlation functions turns out to be 
exponentials: 



Tcit) = (cos iy9(f) cos iy9(0)) — (C0S(^)^ — AcCXp 

Ts{i) = (sin(/3(t) siniy9(0)) = exp 



(34) 



It is important to note that the amplitudes of these correlations (which are equal to the ones of the original 
system as the variables (f are not rescaled) are analytically computable at each temperature and are given 

by 



A, = (cos2 (p) - (cos (p)' = 1 - f c(f ) - c(r)2 , 
As = (sin^ ^) = f C(f ) . 



(35) 
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FIG. 6: Relaxation times of the correlation functions of Eq. 12811 . The symbols are the result of the numerical 
calculation, while the dashed lines are the fits and the dot (dot-dashed) lines is the low (high) temperature expansions. 



where C(T) is the magnetization expressed as a function of T and determined self-consistently by Eq. ((SJ. 
The relaxation times are related to the real ones by Tc^s — jT^^k^T^cs- The limits in which Eq. 134(1 are 
analytically obtained are: 

• High T limit: if T — > oo one can neglect in Eq. ((27|l the term proportional to simp. In this case 
the dynamics is close to a free Brownian motion and the correlation functions are exponentials with 
fc ~ fs ~ 1/f and As ^ Ac^ 1/2. 

• Low T limit: the low temperature limit is obtained by considering ip{t) 7r/2 and i) by expanding 
sm (y9 ~ (/3 m Eq. ^ and ii) by approximating Tsii) - {ip{i)ip{0)) and ~ \[{ip'^{i)(p'^{0)) - 

In the approximation i) the equation of motion l|28|) is easily solved, the correlation functions are 
exponentials and one obtains 

-II ' -I 

" " T ' " 2 ' (36) 

As=f , T, = 1 . 

The expressions for the amplitudes are consistent with Eq. H35|) observing that in this limit C ~ 

-I _ T _ 
^ 1 8 ■ 

The complete T dependence is obtained by solving Eq. ((28|l numerically: the numerical solution has been 
performed using i) the true dynamical system (defined through Eq. (^| and Q), and ii) the effective one 
defined by Eq. H28() . In both cases the reduced variables (t and f\) have been used. The integration of the true 
dynamical system allowed to derive numerically the self and the collective correlation functions (in each case 
for both the sin(/5 and cosiy9 variables). The self correlation functions, that for check have been compared 
to those obtained trought the integration of the effective dynamical systems, were fitted to an exponential 
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FIG. 7: The parameters of the self correlations for k — 2. The symbols are the result of the numerical calculation. 
The full lines are obtained from analytic computation, while the dashed lines are obtained from the polinomial fit 
fEg. I37II on the reduced system and the substitution of T with its value for k = 2 (see text). The agreement is not 
perfect around Tc due to finite size effects. 



decay to derive the parameters Ac, As, and Tg. The simulated systems is composed of TV = 1000 degrees 
of freedom (in the case of the effective dynamical systems, the A'^ independent degrees of freedom have been 
used to collect statistical average over the initial conditions). The equation of motion have been integrated 
by a simple constant stepsize {dt = 0.001) Runge-Kutta method, and the RAND 39] fortran routine has 
been used to generate the gaussian noise fj. At each temperature we performed an equihbration of 5 • 10^ 
integration steps followed by 5 • 10^ steps of data collection. The time history of the variables (f have been 
stored and a multi-step circular buffer scheme has been employed to calculate the appropriate correlation 
functions. Wc found that for any T the correlations are well (but not exactly) described by Eq. 1)34(1 : in Fig.|Sl 
we report some correlations calculated numerically for intermediate T values, together with an exponential 
fit. The relaxation times obtained numerically are reported in Fig. |B1 They have been fitted for simplicity 
with a polinomial function: 

l + PiT + PaT^ + PgTS 

where 'r(O) is given by Eq. ((36|l . The previous expression reproduces the correct high and low T limits. The 
values of the parameters Pi are: 

fc : Pi = Pa = 0.90 Pg = 6.28 , , 

(38) 

fs-. Pi = 1 P2 = 1-66 P3 = 6.28 ^ ' 

The relaxation times obtained by the numerical calculations are reported in Fig. |^ together with the corre- 
sponding fit and the high and low temperature expansions. Having an (exact) expression for the amplitudes 
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FIG. 8: The parameters of the self correlations for fc = 3 (with the same notations used in the previous figure). The 
vertical lines marks the transition temperature To; the numerical calculations can be done also in the metastable 
phases around Tq. The white dots mark the point at which the two solutions with C disappear. 



(Eq. H35|l ) and another one for the relaxation times (Eq. (|37() . extrapolated from numerical data) we can 
discuss the behavior of the correlation functions for any value of k by substituting in these expressions 
f - ^ 



2. k = l 

For k = 1 we have T — AT and Tc^s = ■Jt^c.sj so that the behavior of the relaxation times is obviously the 
same as in Fig. The amplitudes, that we do not report, are analytic functions of the temperature. 



3. k = 2 



As reported before, for fc = 2 a second order phase transition takes place at Tc = A. The parameters A and 
T for this case are reported in Fig. \7\ We have that T ^ oo by approaching the phase transition from the 
magnetic phase, so Acs — > ^ and t^^s ~^ l/Tc — 1 / A. at the transition. In the paramagnetic phase one has 
C = 0, so that Eq. H27|) reduces to a free Browian motion; then Ac = As = 1/2 and Tc = Ts = 7/r exactly 
for all temperatures above Tc- We obtain then that A and r are continuous functions of temperature but 
their derivatives have a discontinuity at Tc. 
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4. fc > 3 

For > 3 a first order phase transition takes place. The parameters A and r are reported in Fig. |H1 In this 
case T does not diverge approaching the transition temperature, while the behavior in the high temperature 
paramagnetic phase is the same as for k — 2. Then, obviously, A and t are discontinuous at the transition. 



VI. COLLECTIVE DYNAMICS 



To study the collective behavior of the system we introduce the correlation function of the magnetization z. 
It is defined by 



•^(^'^'^ = ^Yl [(e*[^'(*)-^^-(*')l) - (e'^»(*))(e-*^^- (*'))] = N [{z{t)z*{t')) - {z{t)){z*{t'))] 



(39) 



Again, with the aim to obtain simple exponential behavior of the correlation functions, it is convenient to 
separate the contributions coming from the real and immaginary part of the magnetization, defining 



Gc{t,t') ^ ^ ^ [{cos (f ^{t) cos ip jit')) - {cos (pi{t)) {cos (fjit'))] , 
Gs{t,t') = ^^(sin^,(t)sin^,(0) . 



(40) 



These correlations are related to the Gaussian correction around the thermodynamic limit, i.e. the leading 
correction in 1/A^ for — > oo. In appendix El we derive a general expression (Eq. !KJT^ ) that relates G{t) 
to the self correlation function F{t) defined in the previous section. Using this expression and assuming 
that i) the magnetization is always real and ii) Fc^s{t) are given by Eq. (|34|) . Gc,s{t) turn out to be also 
exponentials and are given by: 



Gc{t) = Zc Ac Clip 



Gs{t) = Z^^exp 



where 



Zr = 



1 



1 - /3Afc(fc - l)C(/3)'=-2^ 



Zs = 



1 



1 + /?Afc(fc - l)C(/3)'=-M 



As expected in absence of interactions, for fc = 1 we have Zc = Zg = 1, so that G{t) 
temperatures. 

In order to treat the k > 2 case we define 



Gc,s(0 =^?aexp 



t 



with 



A,. „ — Zr.sA, 



c.s-^c.s , 



(41) 

(42) 
F{t) at all 

(43) 
(44) 



From the analytic expression for A^^si given by Eq. (|35|l . we get an analytic expression for the constants 
Zc,s- Exact expression for the parameters A'^^, that are connected by the fluctuation-dissipation theorem 
to the susceptibilities of the magnetization, can then be derived. The relaxation times of G{t) are obtained 
from the one of F{t), that we studied in the previous section, using Eq. H37() . Then, their expression is not 
exact but derives from the numerical data on the reduced system that we defined in the previous section. 
To emphasize this, in Fig. |51and^|we report as a full line the exact expressions and with a dashed line the 
expressions derived using the numerical solution of the reduced system and Eq. H44|l . 
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FIG. 9: The parameters of the collective correlations for k = 2. As in the previous figures, the symbols are from 
numerical computation, the full lines are obtained analytically, while the dashed lines are obtained using the fitted 
expression (Eq. 13711 ) in Eq. 1441 . In the inset the divergence at Tc of the susceptibility and of the relaxation time 
related to the real part of the magnetization is evidenced. 



A. fe = 2 



The parameters A'^ and r*^ given by Eq. (|44|l for fc = 2 are reported in Fig. El From the symmetry structure 
of the model, wc know that it can acquire a spontaneous magnetization in the directions — and — n, 
but not in the orthogonal direction. Then we expect a divergence in the amplitude and in the relaxation 
time of Gc{t) but not in the same parameters of Gs(t). In the high temperature paramagnetic phase we 
have Ac = As = 1/2 and Tc = Ts = 7/r (see the discussion of the k = 2 case in the previous section); 
then Zc = 1/(1 — (3 A) and Zg — 1/(1 + /3A). From these expressions and Eq. (|44|l we get an expression 
for A^ and r*^ in the paramagnetic phase (full lines in Fig. EJ. It is easy to see that close to Tc one has 
Af ~ \T—Tc\~^ and ~ \T—Tc\~^. The same behavior is obtained approaching the transition temperature 
from below, as one can easily check remembering that for T T~ one has T — > 00, Acs 1/2, Tc^s l/T. 
We obtain again the classical mean field critical exponents for the universality class of our model. In the 
low temperature phase, as previously discussed, the expression for the relaxation times in not exact, and is 
reported as a dashed hne in Fig. El 



B. fc > 3 



The parameters A'~^ and r*^ given by Eq. (|44|l for fc = 3 are reported in Fig. ^1 As previously shown, for 
fc > 3 the model undergoes a first order transition at Tq. The low temperature phase is metastable up to a 
certain temperature Tc (see Fig. ^ where it disappears, while the high temperature phase is metastable down 
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FIG. 10: The parameters of the collective correlations for fc = 3 (with the same notations used in the previous figure). 
In this case and diverge in the metastable phase when the metastable solution disappears (white dot in the 
figure). 



to T = 0. In the high temperature phase we have, as in the k — 2 case, Ac = Ag = 1/2 and Tc = = 7/T. 
But for fc > 3 we have Zc = Zg — 1, so that G{t) = F{t). We note that A'^^ does not diverge in this phase, 
and the same happens to t^^ that diverge only for T ^ 0. In the low temperature phase it is easy to see 
(substituting Eq. in Eq. that Z^ — 1/k; this happens for all fc if C 7^ 0. Then, Af and are 

simply proportional to As and respectively, and do not show any anomaly. The behavior of Zc is more 
interesting: using its definition given in Eq. 1)42(1 . Eq. JHJ and Eq. 1)35(1 . one can show that Zc ^ 00 when 
T Tc ■ Thus, A'^ and diverge in the metastable region of the paramagnetic phase approaching the 
temperature at which the paramagnetic phase itself disappears. Again, the expression for t*^ is obtained 
using Eq. , and is reported as a dashed line in Fig. 



VII. CLOSEST SADDLES TO EQUILIBRIUM CONFIGURATIONS 

The important role that stationary points (saddles) of the Potential Ene rgy Surface (PES) plays in the 
dynamics of various systems has been clarified recently [211 123. 123L I251 1261 [21 . 28, 34] . From previous 
works it seems that, in order to describe the equilibrium dynamics at a given temperature T, it is sufficient 
to know the properties of some of them, that have often been called "closest saddles to the equilibrium 
configurations at temperature T" To locate these particular stationary points, two main strategies have 
been adopted: 1 ) defining in a proper way a "distance" in phase space and, given an equilibrium configuration, 
looking at the stationary point that has minimum distance from this configuration; 2) partitioning the phase 
space in "basins of attraction" of stationary points via an appropriate function that has a local minimum 
on each stationary point. While the first approach has been exploited analytically on some disordered spin 
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models | ^4l. t he second one has been extensively used in numerical simulations of simple model liquids 
I2II2II2II2II2II2II23 

The problem with definition 2) is that one has to define a function W such that each stationary point of H 
is a local minimum ofW and each local minimum of W is a stationary point of H. Otherwise, while looking 
for the closest saddle starting from a reference configuration one can remain trapped in some local minimum 
of W that is not a stationary point of H. It has been shown in 24] that this possibility effectively arises in 
the majority of the cases if one chooses, as usual in simple hquids, W — | ViJp. 

To compare the two methods avoiding the difficulties of the numerical computations, we tried to find some 
models in which the minimization of a function W with the desired property could be analytically performed. 
In the fc-trigonometric model this function can correctly be chosen as = |Vi7p, as one can check directly: 
in fact, all the minima of W correspond to stationary points of H. In this section we present a general 
method for the minimization of W = |V-ff p, that can be probably extended to treat a large class of mean 
field systems without quenched disorder. We apply this technique to our model and we show that definitions 
1 ) and 2) give in this case exactly the same result. Note that the idea on which the method is constructed can 
be used also with W functions different from the one chosen here, even if the practical calculation might be 
difficult depending on the particular form chosen for W. Future work will be devoted to apply our method, 
if possible, to disordered systems like the p spin disordered model studied in [33 |. 



A. Definition of the relevant quantities 

In section El we studied the properties of the stationary points of the PES that are independent from the 
statistical measure that describes the system at temperature T (i.e. they are independent of temperature). 
From this study we were able to find a relation between the energy of a saddle and its order and to compute 
the number of saddles, expA^CT(e), located at a given energy e. Now we want to calculate the energy of 
the "closest saddles to equilibrium configurations at temperature T" . The procedure used to calculate this 
quantity is the same used numerically in pit : we consider an initial configuration extracted from the 
Gibbs distribution at temperature T, and we perform a minimization of 



W = \VH\' = [(z z*r-' - ^iz'"-' z,)] , -2 = ^ E e^'^' . (45) 



that leads to the "closest saddle to the initial configuration". Finally, we average over the equilibrium 
distribution of initial data. 

The minimization of W is performed using the dynamical system 

dW 

1<P^=-J^ , (46) 

that is completely analogous to (QJ with H-^W and T=0. We want to calculate the energy of the configu- 
rations in the limit of infinite time starting from a Gibbs ensemble, i.e. 

es{T)^^YiTa{H{t))w ■ (47) 

From this quantity, we obtain the order of the saddles as a function of temperature using Eq. (|20|l and the 
"configurational entropy" of the saddles that is given by cr(T) = a{es{T)). 
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B. Effective dynamical system 

Using the same argument presented in appendix ^ for the real dynamics, it is possible to show that the 
dynamical system (|46|) is equivalent to the single-particle one given by 

at) = (e'^(*)) , 
U2(t) = (e2^'^(*)) . 

... 

The derivation is sketched in appendix^ In this case we will calculate the averages over the distribution of 
initial data setting the noise to zero. As the system is mean field, the correlations between different degrees 
of freedom vanish in the thermodynamic limit and the Gibbs distribution can be written in the form 

piM) = n ^(^^) = n = n — ? — ' (49) 

% % % 

where C(/^) — — 0) is as usual the equihbrium average magnetization. The problem is then reduced to 
the calculation of 

e,(r)= lim(H(</.))w , (50) 

t — >oo 

using the dynamical system and averaging over the distribution of initial data. 



C. Calculation of the energy of the closest saddles 



Again we assume that the system is not magnetized or that the mean magnetization has zero phase. So, we 
assume that and C2(i) are real functions of time. In this case the system (|48|) becomes 



>V(</j,i)- A2fc2 (^{tfk-z [ (I-C2W) cos(^-iCWcos2^] 



I C(i) = (cos(p(t)) 
[C2W = (cos2^(t)) 



(51) 



We want now to show that the first term in W can be neglected, at least in some limits. For fc = 1 the 

A 2 

first term disappears and one simply has W = — ^cos2(/3. For /c > 2, in the paramagnetic phase one 
has C(0) = C2(0) = 0, then W = and the closest saddle is the starting configuration itself. In the low 
temperature phase we know that at the initial time 1 — C2(0) = 2(sin^ ip) = fcAC(o)'="^ • Then we can neglect 
the first term with respect to the second one (at t = 0) if 

2T 1 T A" — 1 

where the 4 in the right side comes from the fact that the derivative of the second term is proportional to 
2sin2(/5 ~ 4(^ while the derivative of the first term is proportional to sin(^ ~ ip. This condition is clearly 
satisfied for low enough temperature because C,{T = 0) = 1. It is easy to check (see Fig. that for fc = 3 
the inequality is satisfied up to the transition temperature Tq. Obviously for fc = 2 it cannot be satisfied 
close to Tc where C '^-^ 0. 

If one can neglect the first term at t = 0, it can be neglected at all subsequent times, because during the 
minimization of W both C(t) and ^2(0 increase, and C2 ^ 1 rapidly. To give an argument, let us neglect again 
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FIG. 11: Saddle energy Cs versus the thermodynamic energy e for k — 1,2, 3. The dashed line is the unstable region 
that corresponds to the solution of Eq. |H|that is a maximum of the free energy. 



the first term in W; then the minima of the pseudopotential are such that cos 2ip ^ 1, so that increasing time 
C2{t) will move toward 1 that is its infinite time limit. C will increase due to the fact that the final energy 
is lower than the initial one, as we will show below. Surprisingly, while we expect all these approximations 
to work only at low enough temperature, they give the correct result in the whole temperature range, as we 
checked numerically. 

In the approximations discussed before, the system H51() becomes of the form: 

r <^ = _ VW - -A^PCity''-^ sin 2ip , , . 

\CW = (cos^W), ^^^^ 

We are interested in the infinite time solutions of Eq. (|53|) . ip{t—>oo\ipo), as a function of the initial conditions 
(fio = (p{t = 0). Without explicitly solving Eq. (|53|l . we observe that the sign of VW at fixed (p does not 
change during time, due to the fact that the time dependent factor in Eq. 1)53(1 is always positive. This 
implies that the specific time dependence of C does not affect the final point ip{t oo\ipQ) reached from a 
given initial condition, rather it controls the rapidity of approaching this final point. It it easy to see that 
the solutions are: 

r if ^o^(f,f) , 

(p{t->oo\ipo) ^ < (54) 
[ TT if ^oe (f,f ) . 

The solution = is a minimum of W corresponding to a minimum of the effective potential energy 
(Eq. 0), while the solution </? = tt is a minimum of W that corresponds to a maximum of the effective 
potential. The energy of the saddle is then obtained as es{T) = A(l — C(oo)'^), where C{oo) is the average 
of cos(p{t~^(X)\(pQ) over the distribution 149|) of initial data: 

C(oo) = J dipo Viipo) cosip{t->oo\(po) = J dtpo V{ipo) sgn(cos(^o) = /"(/gAfcC'^-i) ' ^^^^ 

where we have introduced the modified Struve function of order 0: Lo(ct) — 2tt^^ J^^^^ dip smh{a cos ip). 
In Fig. 12 the saddle energies eg (dashed lines) are reported as a function of temperature for k = 1,2,3: 
qualitatively reproduces the shape of the potential energy e, and it is always below e, but coincides with e 
in the paramagnetic region for k>2, as previously discussed. The map vs. e is shown in Fig. ^2 where one 
observes that, when the parameter k increases, the energy difference between instantaneous configurations 
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FIG. 12: Saddle order as a lunction ol temperature for k — 1,2, 3. 



and saddles becomes more and more pronounced. From Eq. H20|l we obtain the saddle order n as a function 
of temperature: 

n(T) = i 

This function is reported in Fig. E|for fc = 1, 2, 3. At low temperature (high /?) Eq. H56|) is approximated by 

which corresponds to an Arrhenius behavior. We have shown in j4ll | that the behavior of ?t.(T) is related 
to that of the diffusion constant D{T), as previously found numerically in simple model liquids j2ll |. and 
that the energy barrier appearing in the Arrhenius low temperature expansion of D{T) is exactly the energy 
difference between saddles of order 1 and the underlying minima. 



1 - 



£o(/3AfcC"-i) 



(56) 



D. Distance of the closest saddle to the reference configuration 

In this section we show that the two possible definitions of "closest saddles to equilibrium configurations" 
that we discussed above are coincide in our model. To this end, we apply the method introduced in [34] to 
our model. We compute the quantity 

I r „-f3H{<p) r 

a{T;es,d) ^-J d^,-^^\og J d^j, S{H{-^) - Ne^) S{d,H{i:)) det H(^) <5 (d - d((^, ^)) , (58) 

where ~ didjH is the Hessian matrix and d{ip,'il)) is some distance function between the two config- 
urations ifi and ipi. The argument of the logarithm is the number of stationary points of energy eg and 
distance d from the reference configuration ip (see appendix A or Ref. [3^ for a detailed discussion). Then 
the logarithm of this number (divided by N) is averaged over the equilibrium distribution at temperature T 
of the reference configuration. 

Using this quantity we can provide a definition of "closest saddles to equilibrium configurations" [3^ : in fact, 
let the temperature be fixed (and neglect the explicit dependence on it of all the quantities) and consider 
(T(es, d) as a function of d at fixed e^. This is the number of saddles of energy Cs and average distance d from 
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FIG. 13: Mean overlap between instantaneous configuration and the closest saddles (see text) as a function of 
temperature for k — 1, 2, 3. 



the typical configurations at temperature T. We expect that for too small distances this quantity will be 
zero, because there will be no saddles of energy Cs at too small distance from the equilibrium configurations. 
So we can define d{es) as the value of d at which cr(es,d) goes to zero: a{es,d{es)) = 0. Then d{es) is the 
minimum distance at which one can find saddles of energy e^. Now we can minimize c?(es) with respect to 
Cs'. the value Cs of Cs such that d(es) is minimum will be the energy of the closest saddles to the equilibrium 
configurations, while c?(es) will be the average distance from these saddles and the equilibrium configuration 
themselves. 

In our model the distance function can be defined as 



(59) 



In fact cos((pi — ipi) can be interpreted as the scalar product of the unitary spins represented by the angles 
Lpi and ■(/'i, so that q{ip,4') is the overlap between the two configurations. The calculations are reported 
in appendix ^ the result is that the energy of the saddles is given by the same expression derived in the 
previous section: 



£o(/3AfcC^-i) 



where ^ is given by the thermodynamics, see Eq. |(SJ), while the mean overlap is given by 



2tt 



dip Viifi) I cos(^| = (I cos(p\)n ^ 



(60) 



(61) 



where V has been defined in Eq. H49|l and Li{a) = L'^ioi)- Note that this result could also be derived using 
the method of the previous section: in fact, we would have, in analogy to Eq. H5U|I : 



qs{T) = lim (cos((^(i) - <y9o))w = (cos((^(i oo|(^o) - (Po))n(vo) = (I cos((^)|)-h = q{T) 



(62) 



Then, we can conclude that in our model the minimization of W starting from an instantaneous configuration 
equilibrated at temperature T leads, on average, to that stationary point of the PES which has maximum 
overlap with the starting configuration itself. Moreover, we can calculate the average overlap (or distance) 
between the equilibrium configurations and the associated stationary points, that is reported in Fig.^las a 
function of temperature for k — 1,2,3. 
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VIII. CONCLUSIONS 

In this paper we presented a study of the thermodynamics and the dynamics of a very simple mean field model 
of N variables interacting trough a fully-connected /c-body trigonometric term. In spite of its simplicity this 
model shows an interesting behavior undergoing second or first order phase transitions depending on the 
value of k. All the results that we present here have been obtained analytically (except for the computation 
of the relaxation times of the self correlation function that has been done numerically). In particular, 
the dynamics is analyzed in detail, and we find a relation between the single-particle dynamics and the 
collective one in terms of a Dyson-like equation that links the self and collective correlation functions. Using 
this relation we studied, among other dynamical properties, the critical slowing down around the second 
order phase transition or close to the stability limit in the case of the first order phase transition. A 
rather complete characterization of the geometry of the Potential Energy Surface (PES) has been achieved. 
This allowed us to relate some of the macroscopic properties of the model to the PES characteristics: the 
thermodynamic singularities (phase transitions) are located at the same energy of the geometrical ones, and, 
in the dynamics, the low-temperature behavior of the diffusion constant is similar to the behavior of the 
order of the saddles visited during the equilibrium dynamics. The concept of "saddles visited during the 
equilibrium dynamics" (or "generalized inherent structures") has been widely used in the literature, even if 
an unambiguous definition of them has not yet been found. In this respect, we compared two definitions 
that have been used in the past, and we found that they give exaclty the same result when applied to the 
fcTM. This result supports the use of PES to analyze the behavior of interacting systems and suggests that 
the present analysis could be applied to other interesting systems. 
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APPENDIX A: TOPOLOGICAL PROPERTIES OF THE SADDLES IN MEAN FIELD MODELS 



In this appendix we will calculate the number of stationary points of energy e = E /N and their order u sing 
a general method that works well for all mean field models. Although this method is well known |3 l42l | 
and in our model the same results can be obtained in a simpler way, it is useful to recall it here because it 
will be generalized in appendix El to calculate the number of saddles of energy e located at a given distance 
from a reference configuration equilibrated at a given temperature T . 
We introduce the quantity 

N 

X{E) = / d^, &{H - E) S{d,H) detH - ^(-l)^Ar,(£;) (Al) 

i/=0 

where H is the Hessian matrix (Hy = didjH), Mu{E) is the number of stationary points of H of order v 
and energy E, and a product over the index i is omitted. The last equality is easily checked by noting that 
the function detH 5{diH) is equal to (—1)'^ if integrated in a small volume around a stationary point of 
order v. Using the relations 



do.H^ .i:..^^". 



Vj 



(A2) 



where {77, 77} are Grassman variables [36l | we have 

X{E) = J ^e'^^ J dif,dfj,drj,dX, exp [-PH + tj, d^d.H % + 2mXAH] (A3) 
Introducing the superfield 

(j)^(e, 0) = ip, + me + fi,e + 2mxM (A4) 

where 9 and are two other Grassman variables we have 



X{E) = j ^e^^ I Vcb. exp 



ded9{l - l3d9)H{(l)) 



(A5) 



The last equality is easily checked remembering that 

H{(j)) = H{p) + (?7,0 + nfi + iX,e§)d,H{ip) + irj,f],ee)d,d,Hiip) (A6) 



due to the fact that 9'^ — 6^ ~ 0. In mean field models we can evaluate the integral (|A5|I at the saddle-point; 
form Eq. IjAip we see that the integral will be dominated by a particular value of ly: 

X{E) ^M^o. {-\f^''^M,^E){E) ^ e^-(^) (A7) 

so that 

lim 1 log x{E) = lim 1 ( log A^.f^;) [E] + i^E)) (A8) 

Then we expect that at the saddle point the real part of a{e) will be the logarithm of the number of saddles 
located at energy e, while its imaginary part will be the order of these saddles [s^. We will now calculate 
explicitly all this quantities in our model. 

The hamiltonian of the fc-trigonometric model is written in term of the the variable z = e^'^^ in 

Eq. (QJ. This variable has a real and an imaginary part. As we want to evaluate Eq. ljA5p at the saddle 
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point in z, we will need to consider the real and imaginary part of z as complex variables themselves. To 
avoid confusion, it is convenient to use another imaginary unit I and z = e^*'' = SRz + I'^z. When 

we will consider and 3z as complex numbers themselves, we will use the notation Ikz = He^z + ihaStz. 
Setting A = 1 we have (neglecting all the constant prefactors) 



X 



X(e) = j dpe^^ jvzd ^Nz - ^ e^"^' j J V(i)i exp j dMe{l - I366)N{1 - ^z'' 

= J (1(36^1^" j PzPiexp I j Sde [(1 - /36»^)(1 - ^z'') + ^{z iz)] + log A'(ji) 
{iz) = jv^ exp - J dm SR(e^"^ iz 



(A9) 



where we introduced the supervariables z = zo + Zi9 + Z2O + zzOO and z and the superdelta-function 



5{z) = Jt^z exp j dOde ^{z iz) 



(AlO) 



In the definition of 6{z) and of X{iz) wo note that both imaginary units appear: / serves to select only the 
component cos — 9£ sin c6 in the product e^'^z, but 3£ and (f) are themselves complex superfunctions 
with respect to i. We will also rotate the integration path on and 9i in the complex (i) plane (that is 
equivalent to the substitution iz ^ z); this is irrelevant because at the end we will look for the saddle point 
in the whole complex (i) plane. We obtain then 



Xie) 



j d(3T>zT>z e^'^^^'^''^ I "'^ 



a{z, z,(3\e)=Pe + j dOdO [(1 - /36»6')(1 - ^z^) + ^zz] + \ogX{z) 



(All) 



To solve the saddle-point equations we will assume that at the saddle-point: i) all the fermionic components 
vanish (zi = Z2 = z\ = Z2 = Q) and ii) the 9 part of the bosonic components is always {^Zo = ^z^ = ^zq = 
= 0). The first assumption is standard in this kind of computations. The second one is a consequence 
of the symmetry structure of the model: one can always choose the magnetization z = Tli^^'^^ such 
that its imaginary (in the / plane) component is zero. Firstly, we will evaluate X{z) in the case in which z 
has the form that we have assumed above. We get 



X{z) 



f 

Jo 



dip e~'^3cosv S{zosmip)zocosip 



V(j) e 
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If io is real and different from we get 

Then we obtain 

(t{z, z,(3 I e) = Pe + zoia - 
The saddle-point equations are then 



^oizs) = -2 sgnzo sinhz3 
f Z3Z0 - /3(1 - zo) - kz^'^^zs 



log(-2 sgnzo sinhia) 



(A12) 

(A13) 
(A14) 



e = 
zo 
Z3 
Zo 



1 



kz, 



k-l 



-(3kz. 



fc-i 




(A15) 



tanh Z3 



25 



Substituting in a wc obtain 

a{e) = Z0Z3 + log 2 sinh % + log sgnzo (A16) 

Now, if k is odd or if k is even and e < 1, zq is real and given by zq — {1 — e)^/'''. Then zq is real and 
Eq. (|A13|I is correct. Note also that zq is positive (because fc — 1 is even or e < 1) so that the last term in a 
is 0. Recalling that atanhx = ^ log one has 

, f I \ I , I - Zf) ITT / A 1 

Z3 = atanh = o + T 

\ zqJ 2 1 + ^0 2 

It is interesting to note that as £3 is complex while zq is real, form the third of the saddle point equations 
one obtains that f3 is complex at the saddle point; this is a consequence of the srongly oscillating behavior 
of x(e)- Using the relation sinh a; = . ^ we get sinhza = , ' and introducing the variable 

nie) = i(l - zo{e)) = 1 [l - (1 - e)'/'^] (A18) 

we finally obtain 

a{n) — —n log n— (1 — n) log(l — n) — inn (A-19) 

and 

X{n{e)) - (-l)Nn{e)^Nain(e)) (^20) 

This result is consistent with the one obtained in section HVI and with the discussion at the beginning of this 
appendix. The case in which k is even and e > 1 is a little more involved and we will not discuss it here. 



APPENDIX B: SINGLE-PARTICLE DYNAMICS 



In this appendix we will show that one can use the effective dynamical system H25f) to compute one particle 
quantities (see Eq. (|22|l l. We will restrict to the observable e*'''^*) and its n-times correlations; the other 
observables are linear combinations of this one (via a Fourier expansion). We will use the formalism of the 
generating functional in its supersymmetric formulation as presented in |16t ..17. ..36.J , and a notation similar 
to the one of appendix A. 



1. The generating functional 

The generating functional of the correlation functions can be written as [l^ IT^ Is^ 

Z[hit)] = J V(j)i exp ^YlJ '^^ '^*(") ^^^^ '^»(") ^ ^^^^ da h{a) e'*^ 



(Bl) 



where 6*, are Grassman variables, (j)i{6,6,t) is a time-dependent superfield, see Eq. ljA4|l . da — d6d6dt, 
h{a) = h{t)ee, h{t) = m{t) + /9/i(i) and 



d9de ' dOdt ' dOdedt 
In fact it is easy to check that the self-correlation functions (|32|l can be written as 

^, r S^Z 5Z SZ 6Z 5Z 

r (t,t ) 



sm{t)sm{t') s'^h{t)5'^h{t') 6m{t)sm{t') s'^h{t) 5'^h{t') 



(B2) 



(B3) 



h=0 



using the symmetry under permutations of the ipi. Defining formally the operators 

6 6 5 6 (5 



Shit) 6m{t) 53/i(t) dh*it) \Sh{t) 



one can see that 



F{tX) 



sz sz 



6h(t)5h*(t') 5h{t)5h*{t') 
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(B4) 



(B5) 



h=0 



using {sh\{Lpi{t)) cos{ipi{t'))) = and {sm{ipi{t))) = because of the symmetry ip — s- —ip, as we have already 
discussed before Eq. 132|l . We can generalize this relation defining the supercorrelator 



F{a,b) ^ 1^ [(e^(^>(«)-^.W)) _ (e^0.(-))(e-^^-('')) 

i 

s^z sz sz 



5h{a)Sh*{h) Sh{a)Sh*{b) 



(B6) 



h=0 



2. The saddle-point equations 

Substituting Eq. (QJ in Eq. (|Bip and introducing the supervariable 

i 

we get 



(B7) 



Z[0] = 1 = J Vcjj.Vz S ^Nz{a) - ^ e'-^'^") j exp ^ T{(j),) - NA J da {l~ 3?z(a)'=) 

where T(0) — ^ J da (pA'-^^'cj). Using the integral representation of the (5-function 

S{z) — J Pzexp J da 5ft(z(a) iz[a)) 
we get, rotating as usual the integration path in the z plane, 

Z[0] ^ J V(t)iVzVz exp J da z{a) z{a) - j da z{a) e^'^'^"^ + ^T(0,) - iVASR J da {1 - z{a) 



(B8) 



(B9) 



VzVz expN 



VzVz expNL{z,z) 



/ da z{a) z{a) - A5R / da (1 - z(a)'=) + logZ[z] 



where we defined 



Z[z] = J V(j) exp r(0) - J daz{a) e^-^^") 



(BIO) 



(Bll) 



By comparison with Eq. (|Bip we see that Z[z\ is the generating functional for the dynamics of a single degree 
of freedom with energy 



(B12) 
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We can evaluate the integral in Eq. I|B10|I by a saddle point, as usual in mean field models; we get, using 
again the formal operators l|B4|l 

= z{a) + /^kz{af-^ = =^ C(a) = -AfcC(a)'=-i 

(B13) 



5z*{a) 



6L 



5z*{a) 



z{a) + 



5z*{a) 



log Z[z] = 



C(«) 



'mo 



where we defined C and C. as the saddle point values of z and z respectively. So, in the thermodynamic limit 
the dinamics of the system is equivalent to the one of a single degree of freedom with hamiltonian 



Ida) = (e^^(")) 

Setting 6 = 6 = 0, wc get the effective dynamical system for the variable if. 



(B14) 



(B15) 



Self correlation functions 



We want now to show that it is possible to use the effective hamiltonian, Eq. l)B15p . to calculate the self 
correlation function given by Eq. (|B6|) . In fact we obtain in the same way as we obtained Eq. (|B10|) 



Z[h] 



Vz'Dz exp 



N ^ da z{a) z{a) - NA^ da {I - z{a)'') 



+ (TV - 1) logZ[z] + logZ[z + h] 
DzDz exp 



(B16) 



so that 



F(a,6) 



NL{z, z) + logZ[z + h]- logZ[z] 
log Z[z] 



VzVz 



■ exp [NL{z, z)] 



(B17) 



5z{a)5z*{h) 

In the termodynamic limit we know from the saddle-point equations that the integral is dominated by z = C, 
z = C,, and that 

expiVL(C, ~ Z[Q\ ~ 1 (B18) 

so that 

i!]^l£lfLl ^ (e^W(a)-0(^>))^ 



F{a,b) 



5z{a)6z*{h) 



iHiO 



T(a,b) 



(B19) 



The self correlation function is then equal to the one calculated for a single degree of freedom using the 
effective hamiltonian 1B15II . The same arguments is extended to n-times correlations by differentiating 
n-times. 



APPENDIX C: COLLECTIVE DYNAMICS 



In this appendix we will derive a general relation between the self correlation functions and the collective 
ones. The latter vanish in the thermodynamic limit, and are related to the Gaussian corrections around the 
mean- field saddle point that we studied in appendix IbI 
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1. Gaussian corrections to mean field 



We want now to calculate the correlation function of the "magnetization" (multiplied by N in order to have 
a well defined thermodynamic limit) 



G{a,b) = N[{z{a)z*{b)) - {z{a)) {z* {b))] 



(CI) 



by expanding L at second order around the saddle point; to do this, it is not possible to use the operators 
(jB4ll . We have to separate the real and imaginary part of z and z. We get, calling Az^ = [^{z — C), Q{z — 
C), 5R(2 - C), 3(z - C)] with cr = 1, 2, 3, 4: 



L(z,z)=L(C,C) + 



1 



/ da db Azcr{a) - — , 



)5z^,{b) 



Az^,{b) 



Then, at second order around the saddle point, 



P(Az) c>c exp 



N 



j da db Azo-(a) Jaai{a,b) Azo.'(6) 



where the matrix J(a, 6) is given by 
J{a.b) = 



3?w(a) 5{a - b) -5w(a) 5{a - b) 5{a -b) ' 

-^w{a) 5{a - b) -^w{a) 5{a - b) -6{a^b) 

5{a-b) Tc{a,b) 

-6{a-b) Ts{a,b) . 



where 



J^c(a,6) =(cos(0(a))cos((/)(6)))„(^-) - (cos(0(a)))„(^) (cos(<?;.(fe)))„(^) 

(cos(0(a))sin(0(6)))^(^^ = 
Ts{a,b) =(sin((?!)(a))sin(</>(6)))„(^) - (sin((?!)(a)))^(^^ (sin(0(&)))„(^-^ 
w(a) =fc(fc- l)C(a)'=-2 

are calculated on the effective hamiltonian l)B14(l . We have also 

Jaa'ia,b) = Ja'cr{b,a) 

Then defining 

=N{Az,{a)Az,,{b)) 

one has 



J dbJcrcr'{a,b)Gcr'a"{b,C) = (5o-o-"'5(a- c) 



We are interested in calculating 



G(a,6) =7V[(z(a)z*(6)) - (z(a))(z*(6))] = iV(Azi(a)Azi(6)) + iV(Az2(a)Az2(6)) 
= G'ii(a,6) + G22(a,6) 



(C2) 



(C3) 



(C4) 



(C5) 



(C6) 
(C7) 

(C8) 
(C9) 
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By writing explicitly some of the Eq.s (|C8|I and making some substitutions one gets 

Gii{a,c)^Tc{a,c) + j dh Tc[a,h) [3?w(6)Gii(6, c) - 3u;(6)G2i(6, c)] 
Gi2(a,c) = JdbTcia,b) [3fJu)(6)Gi2(6, c) - 3w(6)G22(^ c)] 
G2i{a,c) = - j dhTs{a,b) [3u;(6)Gii(&, c) + 5Ru.(6)G2i(6, c)] 
G22(a,c) =^s(a,c)- ( dh T,{a,b) [5«;(6)Gi2(6, c) + 5Ru)(6)G22(&, c)] 



(CIO) 



(Cll) 



These equations give the collective correlation functions in term of the self correlations evaluated on the 
effective hamiltonian; but using Eq. IjBlQjl we can replace the effective self correlations with the original 
ones, and obtain a relation between self and collective correlations of the original system. 

2. Solution in the equilibrium case 

Probability conservation and causality imply that |l6l llTj 

G^^. (a, b) = G,,- {t, t') + {9' - 9) {9'R,,, {t, t') + {t, t')) 
J'cAa^b) = Cc,sit,t') + - 9){9"RcAt^t') + 9ncAt,t')) 
where 

R{t,t') = R{t',t) (C12) 

In the t — > oo limit the system reaches equilibrium, so that the correlation function depend only on the 
time difference, the response function are related to the correlations by the fluctuation-dissipation theorem, 
and w does not depend on time and is equal to the equilibrium magnetization. From the solution of the 
thermodynamics we know that the equilibrium magnetization does not depend on 9 and 9, so that 

lim w{a) = w (C13) 

t — >oc 

and 

lim G^^. (a, b) = G,,- {t-t') + [9' - 9) {9'R,,, {t-t')+ 9R^^, {t' - t)) (C14) 

t,t — *oo 
t-t'=r 

If we suppose that the correlations decay fast (exponentially) for t — > oo, the values of ti, in the integrals in 
Eq. (jClOp must stay close to ia, tc- Then in the limit in which ta^ tc go to infinity also tb must go to infinity, 
and w becomes a constant also with respect to the fermionic variables. Then we can rewrite Eq. (jC10|l 

Gil = Tc + ^w Gil -^w Tc® G21 
G12 =^wTc® G12 ~^wTc® G22 

G21 = -3w Ts ® Gn -^wTs® G21 ^ ' 

G22 = - Ts ® G12 - 5Rw Ts ® G22 

where 



(Gi ®G2)(a,c) = dbGi{a,b) G2(6,c) 



(C16) 
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If Gi and G2 are of the form I|C11|1 . d ® G2 has the same form 

(Gi ® G2)(a, c) = (Gi (E)R2 + Ri(E) C2){t, t') + (0' - e){e'{Ri ® i?2)(t, i') + 0{R2 ® Ri){t, t')) (C17) 

where is the time convolution product. Equating (for example) the 9' 9' component of Eq. ljC15|) and 
moving to the frequency domain we finally get 

Rii{uj) = 7^c(^^) + 'Siw 7^cMi?ll(a;) - 7^c(w)i^2l(t^) 

i?12(w) = 7^c(c^)i^l2(w) - nc[u)R22{uj) 

i?2i(w) = 7^,(w)i^ll(w) - 7^,(w)i?2l(w) ^ ^ 

i?22(^^) = ns{uj) - 3w 7^,(u;)i^l2(^^) - 5Rw 7^,(t^)i^22H 

and 

7^c(w) 



Riiiuj) = 



p , ^ ^ ' 

1 + TZsiuj) - i-5R^V^(^) 

Using the fluctuation-dissipation theorem, one can transform this relation in a relation between the correla- 
tion functions. 



3. Other simplifications 

In our model one can always choose the magnetization such that 3w = 0. In this case one has 

7^c(w) 



1 — w TZc{lo) 
Tlsiuj) 



(C20) 



1 + w TZs{uo) 

Moreover, we can assume that the self correlation functions are exponentials, so that 



We get easily 



nu^) = J^^^ (C21) 

1 - liOTc.s 



Rii{^) = , ^^'1' Z, = {1- PwA,)-^ (C22) 

1 — lUjZcTc 



The same relation holds for _R22 with the substitution w ~w. From this relations, using the fluctuation- 
dissipation theorem and moving back to the time domain, one gets Eq.s (I41fl . 

APPENDIX D: EFFECTIVE DYNAMICAL SYSTEM FOR THE MINIMIZATION OF W 



The calculation of the effective potential W used in Eq. (|48|l is carried on the same lines of the one presented 
in section 2 of appendix IBl We start from the generating functional 



Z[0] = j V(j), exp da (j),{a) A^^) 0,(a) - j da 



W{(j>) 



(Dl) 
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where is given by Eq. H45|l . Now we introduce 6 functions for the variables z{a) and 2:2(0), and we get 



Z[0]= J VzVzVz2Vz2 expN di J da{z{a)z{a) + Z2{a)z2{a)) - j da 



N 



+ logZ(z, Z2) 



where now 



Z{z, Z2) — / exp 



Ti 



{4>)-^l da ( z{a) e^^(°) + £2(0) e^^-^^") ) 



We find then that 



The saddle point equation are 



W((/), a) = 3? [z(a) e^"^ + ^2(0) e^^"^] 



(D2) 

(D3) 
(D4) 



,5£ 

bz*(a) 

5L 
5z*{a) 

SL 
Sz*{a) 

SL 



z[a) — — , . , , = 



^(«) + 7^1ogZ[z,Z2] =0 ^ C(a) = (e'*^"^ 



C(a) = A^fc^(fc - 1) [ C(a) (C(a)C*(a))'=-' - C2W C(a) 



W(C,C2) 



1 



2,2/-/„\2k-2 



5z*{a) 

1 

^2 a) - — . = 

^2(a) + ^^logZ[z,Z2] = =^ C2(a) - (e''^^"^)w(c,6) 



(D5) 



C2(a)=--A^fcW 



(5z5(a) 
and finally we obtain 

' W(0, a) = A^k^ik - 1) 3fJ{ [ C(a) (C(a)C* (a))''"' - (2(0) ((a)"""^ ] } - ^A^P s)ff{ ^(a)2fe-2 g2/0 | 
C2(a) = (e2^"*('')) 



that after setting 9 — 9 — reduces to Eq. 



(D6) 



APPENDIX E: CLOSEST SADDLES TO EQUILIBRIUM CONFIGURATIONS 



In this section we will derive the result presented in section FV 11 PI We have to compute the quantity 



a{T;es,q) ^ — j dip 



' Z{T) 



log / d^,, S{H{^) - Nes) 6{d,H{i^)) det H(^) 6 {q - q{^, ^)) , (El) 



where q{(p, ip) — N ^ cos{ipi — tpi). To do that, we need to prove a general relation. Suppose we want to 
calculate at the saddle point a quantity Q of the form 



1 



rf^.^(^logA(^)=lim — 



dipi ——A'^iLo) - 1 = lim log / difii -— yl"((p) , 

^ Z{T) J n^o Nn ^ Z{T) ' 

(E2) 

where we used the relations logx = lim„_,o ^ and lim„^o(/("-) ~ 1) = linin— »o log/('i) if /("■) ^n^o 1- 
Suppose also that the energy depends only on some collective parameter as in mean field models: H{ip) — 
Ne{z{ip)), where Nz{ip) — z{ipi) (in our model, z{Lpi) — exp(/(pi)). Then we have 

-fiNeiz) 



1 f „ 
Q — lim log / dz — — 



-l3Ne{z) 



■ lim — — log / dzdz- 
n^o Nn ' 



■ lim log — ; — r 

n-^o Nn ^ Z{T) 



dip, 6{z-z{ip)) A'\ip) 



dip, ^"(^) 



Z{T) 

dzdz e-/3^(<-)-^-(";-.-)) 



(E3) 



where we defined 



s{n; z, iz)=ziz+^ log j dip, e'^^^^ ^('^'^ 
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(E4) 



Clearly s(0;z,?z) is the entropic contribution to the free energy as a function of z, z that we obtain in the 
calculation of the partition function Z{T), so that 

/(T) = --^ log ZiT) = min[e(z) - rs(0; z, iz)] = e(C) - Ts(0; C, C) = /(O; C, C) , (E5) 

where {C{T), C{T)) is the (T-dependent) thermodynamic minimmn of the free energy (note that at the saddle 
point zz = C). Then we have 



n^o Nn J 

We can now expand z = C + n(^^^'> + o{n^), iz = C, + nQ^^ + o(n^) and 



(E6) 



f{n- z, iz) - /(O; C, C) = ^(0; C, C) + |4 (0; C, C) + 1^ (0; C, C) ^ + o{n^) = 1^(0; C, C) n + o(n2) , 



9z 



diz 



dn 



dn 



because by definition of (C,C) we have §^(0; C, C) = 0, ^(0; Cj C) = 0. We get then the final result: 

g = -/3|^(0;C, = 1^(0; CO • 



(E7) 



(E8) 



We have then to calculate (neglecting the term that vanish on taking the derivative with respect to n): 



s(n;C, 63,9) = ^ log / d^, g-^:. Ccos.,. -Q / d.^'^ 5{H{r)- Ne,) 5{m{r)) detY{{r) 5 {q - q{^,r)) 



N ■ 



(E9) 

where from the thermodynamic calculation C,{T) — —fikC,^^^ and C, is given by given by Eq. (jHJ. Using a 
representation analogous to Eq. we get 



s(n;C,e.,g)=^log / d^, g-^-^"™^'^' 



Y\ j ^e^^"^= j V^l exp j 



1 - /?,M)i?(^''') 



<5 iVg-^cos(^,-^f) 



(ElO) 



We will now: i) substitute the expression i/(^'") = N{\ — SRy^), using y instead of z to avoid confusion with 
the thermodynamic variable C; H) insert some (5-functions for ya and the corresponding integral representation 
with a multiplier jja] Hi) neglect all the product and sum signs related to the index a; iv) use the integral 
representation for the (5-function of q with a multiplier Aq. Then we get an expression that has to be 
maximized with respect to all the parameters to get the saddle point value of s(n; C q): 



s{n; C, es,q) = max 

all par 



Y^Paes+Y, J rf^d0[(l-/3ae^)(l-5R2/a)+%aya]+^Aag + log5(C,ya,Aa) 



SiCyaAa) = / dip V^^e^p 



-C cos - ^ J dOde •^yae^^" -^^-^ ~ 



(Ell) 
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As usual, we will assume that: i) there is symmetry between the replicas (ya = y, etc.); ii) y and y are real; 
in) all the fermionic components are 0. Then we get 



sin;C,es,q) = max n[ P{es - 1 + yo) - ky^ 2/3 + ^02/3 + ysj/o + + log5(C, y, A) 

all par 



3 — C cos 



exp - / dOdO {ijo + ysOO) cos * - A cos{tp - tp) 



(E12) 



Now we have to take the derivative of s with respect to n at n 0. By direct computation 



ds 



ct(C; es,q) = max —(0; C, e^, g) = max ^(e^ - 1 + y^) - ysikyo - yo) + myo + M 

all par On all par 



dip 



-(^ cos ip 



2^/0 (C) 



log / exp I 



d6d9 {yo + y?,96) cos * - A cos{ip - tp) 

(E13) 



The interpretation of this expression is straightforward by comparison with Eq. I|A14|I : in fact if we put 
A = we get exactly Eq. (|A14|I . This correspond to integrating a over so the dependence on the reference 
configuration (and hence on the temperature) disappears and we get the number of saddles of energy Cs. 
When A is different from the last term of the previous expression represents the single-particle version of 
(J. Now we can proceed exactly in the same way as we proceeded after Eq. I)A14(I : we take the derivatives 
with respect to (3 and 7/3. This fixes yo = (1 — eg)^/^ and yo = ^^o"^ ^'^^ equals to zero the first two terms 
of a. As we are looking for saddles of energy < 1, we have then yo > 0, and then the dependence on yo 
in the last term disappears. We get 



(t(C; Cs, g) = max 



/ 



y3yo + Aq + I dip V{(p) log(-2 sinh(y3 + A cos ip)) 



(E14) 



where V{ip) is given by Eq. (|49|) . From the equation yo = (1 — CsY^'' we see that yo is the average of cosip 
on the saddles; then we will change the notation calling yo — Cs- Taking the derivatives of a with respect to 
A and ya we get 



G --y dipViip) [tanhu((^)]-i 
q = — dip V{ip) cosip [tanhu(<y9)]~^ 



(E15) 



where u{ip) = 2/3 + A cos ip. Now if we want that ^ e [0, 1], u{ip) must have an imaginary part; but if we want 
to be real, this imaginary part must be constant and equal to tt/2. We will then assume that ys = y + i§; 
note that this is the correct solution for A = (see Eq. (|A17|l ). We obtain easily the parametric relation for 



G(y,A) 



dp V{p) f{t{p)) 



a{y,X)^ I dpV{p)\\og2cosht{p)-t{ip)f{t{ip)) 



(E16) 



. i-C(y,A) 

ZTT 



where t{p) = y + Acostp and f{t) = ^^l^j^- We see that the imaginary part of a is, as expected, equal to 
i7rn(es), remembering the relation between energy and order of each stationary point. In the following we 
will neglect the imaginary part of a. Now we have to maximize q on the curve a — (see section IVII D|) . To 
do that, we start with a simple argument: as A is the field conjugated to q (the relation between A and q 
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is the same as the relation between magnetization and magnetic field in a ferromagnet) we expect that the 
maximum overlap will be obtained in the | A| — > oo limit. In fact, for A — > ±00, we have t{(p) — > ±sgn(cos (p) 00 
and f{t{(p)) — > ±sgn(cos(^), so that 



lim q{y, X) = T / dip V{(p) \ cos ip\ 



(E17) 



As we want q to be positive, we have to choose A —00. We have then 



lim q{y, A) = dip V{ip) \ cos ip\=q , 



lim C{y,X) ^ dip V[ip>) sgn(cos.^) = Cs , 
lim (j{y, A) = . 



(E18) 



Then we have a consistency check of our assumption, that the point reached in the limit A ^ — 00 belongs 
to the curve cr = 0. 

A more accurate argument can be given in this way: one can look numerically at the curve q{Cs), paramet- 
rically in A at fixed y. The curve looks like the one reported in the figure above, and as A moves from —00 
to cxD the point in the (^5, q) plane moves from (Cs, q) to (— Cs, —q) (black dots in the figure) as predicted by 
Eq. I)E18|I . Now it is easy to show that 



lim q{y,\) = TC 

y^±oa 

lim Cs(y,A) = Tl 

y — >'±oo 

lim cr(y,A)=0 

y — >-±oo 



(E19) 



Then for |?/| ^ 00 and A fixed the whole curves collapse on the white dots in the figure above. But it is also 
easy to show that for y 00 and A = —y{l + 6) ^ —00 the point goes on the upper border of the domain 
reported in the figure above, and moves from the left white dot to the upper black dot as S moves from 
to 00 (the other branch, form the left white dot to the lower black dot, is obtained for 6 going from —2 to 
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— oo). On the whole border of the domain we have cr = from Eq. ljE19|l : then (Cs, q) is exactly the point in 
which q is maximum on the curve cr = 0. 
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